function px=pol(a,x)
% Calcolo del valore del polinomio: 
%
% a(1)x^(n-1) + a(2) x^(n-2) + .... + a(n-1) x + a(n)
% 
% nel punti contenuti nel vettore x
%
n=length(a);
px = a(n) + a(n-1)*x;
xn = x;
for i=n-2:-1:1
   xn = xn.*x;
   px = px + a(i)*xn;
end

   